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Abstract 

The full Dicke model describes a system of N identical two level-atoms coupled to a single-mode 
quantized bosonic field. The model considers rotating and counter-rotating coupling terms between 
the atoms and the bosonic field, with coupling constants g\ and g>2, for each one of the coupling 
terms, respectively. We study finite temperature properties of the model using the path integral 
approach and functional methods. In the thermodynamic limit, N — > oo, the system exhibits 
phase transition from normal to superradiant phase, at some critical values of temperature and 
coupling constants. We distinguish between three particular cases, the first one corresponds to the 
case of rotating wave approximation, which g\ ^ and #2 = 0, the second one corresponds to the 
case of g\ = and gi ^ 0, in these two cases the model has a continuous symmetry. The last one, 
corresponds to the case of g\ ^ and gi ^ 0, which the model has a discrete symmetry. The phase 
transition in each case is related to the spontaneous breaking of its respective symmetry. For each 
one of these three particular cases, we find the asymptotic behaviour of the partition function 
in the thermodynamic limit, and the collective spectrum of the system in the normal and the 
superradiat phase. For the case of rotating wave approximation, and also the case of g% = and 
#2 7^ 0, in the superradiant phase, the collective spectrum has a zero energy value, corresponding 
to the Goldstone mode associated to the continuous symmetry breaking of the model. Our analyse 
and results are valid in the limit of zero temperature, — > oo, in which, the model exhibits a 
quantum phase transition. 
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1 Introduction 



The Dicke model is an interesting spin-boson model because, being a simple model, exhibits 
the superradiance effect [1]. This model describes a system of N identical two level-atoms coupled 
to a single-mode radiation field, simplified according to the rotating wave approximation. In 
this context, the super-radiance is characterized as the coherent spontaneous radiation emission 
with intensity proportional to N 2 . Thermodynamic properties of the Dicke model were studied 
in the thermodynamic limit, N — > oo. It is found that, the model exhibits a second order phase 
transition from normal to superradiant phase at certain critical temperature and sufficiently larger 
value of the coupling constant between the atoms and the field [2] [3]. The influence of the counter- 
rotating term on the thermodynamics of the Dicke model also was studied in the literature, [6] 
[7]. Using different coupling constants between the rotating and the counter- rotating coupling, 
it is calculated the critical temperature and the free energy of the model [4] [5]. We call this 
generalization of full Dicke model. Path integral approach and functional methods were used for 
study spin-boson problems, finding critical temperature, free energy and collective spectrum of 
the models, in the thermodynamic limit, [8] [9]. With this approach, Popov and Fedotov [10] 
[11], rigorously calculated the partition function and collective spectrum for the Dicke model 
in the normal and superradiant phase. Relation between the phase transition and continuous 
symmetry breaking in the Dicke model was pointed out in reference [12]. The full Dicke model 
was studied using the path integral approach [13], here the authors find the asymptotic behaviour 
of the partition function and collective spectrum in the normal phase. Using the same approach, 
thermodynamic properties of some other spin-boson models were also studied [14] [15]. 

In this paper, using the path integral approach and functional methods, we find the asymptotic 
behaviour of the partition function and collective spectrum of the full Dicke model in the thermo- 
dynamic limit, N — > oo, in the normal and super-radiant phase. The full Dicke model exhibits 
phase transition from normal to superradiant phase, at some critical values of temperature and 
coupling constants. In our study we distinguish three particular cases. The first one corresponds 
to the case of rotating wave approximation, g\ ^ and g 2 = 0, in this case the model has a 
continuous symmetry, which is associated to the conservation of the sum of the number excitation 
of the N atoms with the number excitation of the boson field. The second case corresponds to 
the model with g 1 = and g 2 ^ 0, in this case the model also has a continuous symmetry, which 
is associated to the conservation of the difference between the number excitation of the N atoms 
and the number excitation of the boson field. The last one corresponds to the case of g\ ^ and 
92 7^ 0, which the model has a discrete symmetry. The phase transition in each case is related to 
the spontaneous breaking of their respective symmetry. For the case of rotating wave approxima- 
tion, and also for the case of g\ = and g 2 ^ 0, in the superradiant phase, the collective spectrum 
has a zero energy value, corresponding to the Goldstone mode associated to the breaking of their 
respective continuous symmetry. The collective spectrum obtained in this paper is valid for the 
zero temperature limit, corresponding to the case of quantum phase transition. 

Practical realization of the full Dicke model in the laboratory was discussed by Dimer et al. 
[16]. Since the radiation frequency and energy separation between the two levels of the atoms 
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exceed the coupling constant strength by many orders of magnitude the counter-rotating terms 
have a little effect on the dynamics. These authors proposed that in cavities with the N qubits, 
only one mode of quantized field and classical fields (lasers), it is possible to obtain an effective 
Hamiltonian equal to the full Dicke Hamiltonian. It is possible to control the parameters in 
this effective Hamiltonian, and it is possible to operate in the phase transition regime. Other 
authors stressed the importance for quantum information technology of experimental realization 
of generalizations of the Dicke model in cavity quantum electrodynamics [17] [18]. 

Quantum phase transition of the Dicke model, in the thermodynamic limit, is studied by 
diagonalizing the Hamiltonian [19]. For this purpose it is applied the Holstein-Primakoff map, 
which represents the total angular momentum of the N atoms by a single bosonic field. These 
author find the collective spectrum in the normal phase. Similar method was used by Emary and 
Brandes to study the connection between the quantum phase transition and the quantum chaos 
in the Dicke model without using the rotating wave approximation [20]. They find the collective 
spectrum of the model in the normal and superradiant phase, as another quatities properly of 
quantum chaos. The relationship between entanglement and quantum phase transition in the Dicke 
model was also studied [21] [22], the authors find that the atom-field entanglement entropy diverges 
at the critical point of the phase transition. Studies of this relationship between entanglement 
and quantum phase transition for others collective models exist in the literature [23]. 

This paper is organized as follows. In section 2, we introduce the full Dicke Hamiltonian and 
study its symmetries. In section 3, we introduce a map between the spin momentum operators 
of each atom, with bilinear forms of fermionic operators, defining the fermion full Dicke model. 
In section 4, we are able to introduce the path integral approach for the full Dicke model, using 
functional methods we obtain the critical temperature and the asymptotic behaviour of the parti- 
tion function in some particular cases of the model. In section 5, partition function and collective 
spectrum of the model are presented in the normal phase. In section 6, partition function and 
collective spectrum of the model are presented in the superradiant phase. In section 7 we discuss 
our conclusions. In the paper we use ks — c — % — 1. 



2 The full Dicke Hamiltonian and symmetries 

The full Dicke model describes a system of N identical two level-atoms coupled to a single- 
mode quantized bosonic field. The model considers rotating and counter-rotating coupling terms 
between the atoms and the bosonic field in the Hamiltonian, with coupling constants g\ and g^, 
for each one of the coupling terms, respectively. Consequently, the Hamiltonian of the full Dicke 
model can be written as 

Q N N N 

H = T g C W W* b + g (»„+ + + (f«- w + . (1) 



2 



In above equation we define the operators = | (<7^ ±ia^), which the operators cr^, a 2 -^ and 
a ti) = a U) sa ^ s fy the commutation relations [oq\, ofU] = 2 e pqr aL^ with p, q, r — 1, 2, 3. Therefore, 
[<7^, cr^] = (j^ and [c^), cfj)] — ^ z ^ a fj)- The 6 and are the boson annihilation and creation 
operators of mode excitations that satisfy the usual commutation relation rules. 

Let us define three different operators. The first one, the operator N, which is defined by 




The second one, the operator AL, defined by 

1 N 

N- = #b (3) 
z i=i 

Finally, we define the parity operator II by 

U = e i7rN , (4) 

with operator iV defined in Eq. (2). In particular case of g± ^ and g 2 — 0, which corresponds to 
the rotating wave approximation case, it is possible to show that [H, N] = 0. In particular case 
of g>i = and g 2 ^ 0, it is possible to show that [H, AL] = 0. And it is possible to show that 
[if, II] = for arbitrary non-negative values of g\ and g 2 . These commutation relations of the 
Hamiltonian with each operator defined above, correspond to symmetries of the model for each 
case. It is interesting to see that, for the case of g\ ^ and g 2 ^ 0, we only have that [H, IT] = 0, 
it means that the system only has parity symmetry. The operators defined by J p = \ Y^=\ 
with p = 1,2, 3, satisfy the usual angular momentum commutation relations. The Hilbert space 
corresponding to the atoms states can be generated by the basis {\j m)} with j = N/2 and 
m = — j, + — each basis state satisfies J 3 \j m) = m\j m) and J 2 \jm) = j(j + m). 

The Hilbert space, which the photon states are defined, can be generated by the basis with 
their elements satisfying tfb\n) = n\n), in this case, n is the number of photons. Now we are able 
to construct a basis for the total system as a tensor product of the above basis introduced, i.e., 
the set |Jn) ® \jm)}. The symmetries mentioned above, are related with conserved quatities. In 
the case of g\ ^ and g 2 = 0, with [H, N] = 0, the excitation number of the system, n + m, 
is conserved. It means that the temporal evolution of a state given by \n) ® \j m) only evolves 
toward another states \n') <g> \ j m') which n' + mf = n + m. In similar fashion, for the case of g\ = 
and g 2 ^ 0, with [H, A^_] = 0, the difference of excitation numbers, n — m, is conserved. When 
g-i ^ and g 2 ^ 0, which [if, II] = 0, the value e l7T ( n + m ) i s conserved. It means that the temporal 
evolution of a state given by \n) cg> \jm) only evolves toward another states \n') <S> \j m') with 
both, n + m and n' + m! being even orn + m and n' + m! being odd. In all mentioned cases, the 
phase transition is related to the spontaneous breaking of their respective symmetries. In further 
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analysis we shall see that, the symmetry associated to the commutation relation [H, II] = is 
discrete, and the symmetries associated to the commutation relations [if, N] = and [H, N_] = 
are continuous symmetries. In cases of continuous symmetry breaking the Goldstone theorem is 
valid, with the appearing of zero energy value in the phase with the symmetry broken. 



3 The fermion full Dicke model 



Let us define the fermion full Dicke model. For this purpose, let us define the raising and 
lowering Fermi operators a\, (3} and fa, that satisfy the anti-commutator relations aia^+a^ai = 
5ij and + (3j(3i = Sij. In this analysis, we use a representation of the operators cr^, av^ and 

a U) ky the following bilinear combination of Fermi operators, ajcti — (3}(3i, and ft\a.i, the 
correspondence is given by 

ofi) — > a l a i ~ PiPi , ( 5 ) 

a (i) — >■ a JA . ( 6 ) 

and 

%) — >• A- • ( 7 ) 

Using this representation given in Eq. (5), Eq. (6) and Eq. (7) in the full Dicke Hamiltonian 
given by Eq. (1), we define the Hamiltonian of the fermion full Dicke model Hp. So that, we have 

Q N N N 

H F = u tfb+-Y: (oioi - fit fit) + 7= E °iPi + * f #Oi) + 7=E( fota iA + b (5\a) . (8) 

i=l V-'' j=i V-'' j=i 

We are interested in studying thermodynamic properties of the system, therefore we must 
find the partition function Z. It is important to note that Hamiltonians H and Hp are defined 
in different spaces. Each operator af appearing in the Hamiltonian H acts on two-dimensional 
Hilbert space, notwithstanding, Fermi operators a\, a,, 0\ and appearing in the Hamiltonian 
Hp act on four-dimensional Fock space. The following property relates the partition function of 
the full Dicke model with the partition function of the fermion full Dicke model: 

Z = Tr(exp(-f3H)) = i N Tr (exp (-pH F - y N^j . (9) 

In this last relation H is given by the Eq. (1), Hp is given by Eq. (8) and the operator Np is 
defined by 

Np = J2(®la t + Pl(3 t ). (10) 
i=i 
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The traces used in Eq. (9) for each Hamiltonians are carried over their repective spaces. The 
relation given by Eq. (9) let us express the partition function of the full Dicke model Z using the 
fermion full Dicke Hamiltonian given by Eq. (8). 



4 The partition function with path integral approach 

In this section we perform calculations in order to obtain an asymptotic expression for the partition 
function Z of the full Dicke model in the limit of N — > oo. For this purpose we use path integral 
approach and functional methods. Let us define the Euclidean action S of the full Dicke model in 
the following form 

S = ^dT^(r)9 T KT) + E(<(^)^(r)+^(r)9 T A(r))^ - drH F {r) , (11) 

the Hamiltonian Hp is the full Hamiltonian for the full fermion Dicke model, which is given by 

Q N / \ 
H F (r) = u b*(r)b(r) + - £ [a*(r) a< (r) - fiKr)fii(r)) + 

+ j= E(<(t)/^(t)6(t) +a i (r)p;(r)b*(r)) + 

n N / \ 

+ E («*W^(r)6(r) + a*{T)^{r)b*{r)). (12) 

Let us define the formal quotient of the partition function of the full Dicke model and the partition 
function of the free Dicke model. Therefore we are interested in calculating the following quantity 

Z_ = JJM exp [S - f p Ig n(r)rfr) 
Zo J[dr ] ]exp(S -§Jgn(T)dT) , 

the function n(r) is defined by 

N 

n(r) = E0*W «iW + A*(t)A(t)) , (14) 
i=i 

5 = S(b,b* ja,^ , (3, f3<) is the Euclidean action of the full Dicke model given by Eq. (11), So = 
So(b, b*, a, a} ', (3, ft) is the free Euclidean action for the free single bosonic mode and the free atoms, 
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i.e., the expression of the complete action S taking g\ — g 2 — and finally [drj] is the functional 
measure. The functional integrals involved in Eq. (13), are functional integrals with respect to 
the complex functions b*(r) and 6(r) and Fermi fields a*(r), «j(r), (5*(r) and fiiij). Since we are 
using thermal equilibrium boundary conditions, in the imaginary time formalism, the integration 
variables in Eq. (13) obey periodic boundary conditions for the Bose field, i.e., b(/3) = 6(0) and 
anti-periodic boundary conditions for Fermi fields i.e., oti{0) = —0^(0) and Pi(/3) = — Pi(0). 

In section 2, we have analysed the symmetry of the model studying commutation relations 
between the Hamiltonian given by Eq. (1) with some operators defining the symmetry. Now 
we are able to analise the symmetry of the model studying the invariance of the action given 
by Eq. (11) under symmetry transformations. In this way, let us introduce the following field 
transformation 

b(r)^e^b(r), a(r) -+ e ie a(r) , 0(r) e^/3(r) , 

6*(r) ->• e-^b*(r), a*(r) ->• e^ a*(r) , /3*(r) ->• e"** /3*(r) . 

In the case of #i 7^ and g 2 = 0, corresponding to the case of rotating wave approximation, its 
respective action is invariant under tranformation given by Eq. (15), taking 7 = — 0. In the case 
of gi — and g 2 7^ 0, its corresponding action is invariant under tranformation given by Eq. (15), 
taking 7 = <p — 6. Finally, in the case of g\ 7^ and g 2 7^ 0, its corresponding action is invariant 
under tranformation given by Eq. (15), with 7 = $ — = or 7 = 6* — = 7r. In the two first 
cases, the case of g\ 7^ and g 2 = 0, and the case of gi = and g 2 7^ 0, their respective actions 
are invariant under continuous transformation, U(l), of the boson field b(r). In the case of g\ 7^ 
and g 2 7^ 0, its action is invariant under discrete transformations, Z 2 , of the boson field b(r), i. e., 
6(r) -)• 6(r) and 6(r) -)• -6(r). 

Following with the purpose of calculating the quantity given by Eq. (13), let us use the 
following transformation 

«i(r) -)• ew ai{r) , a*(r) -)■ e v l a*(r) , 
ftW-^eft'ftW, P;(T)^e-^P;(r). 

With this last transformation, the term n(r) appearing in Eq. (13) can be dropped. Therefore, 
applying the tranformation given by Eq. (16) into the expression given by Eq. (13), we obtain 
that 

Z_ I[dy}e s 

Z I[d V ]eSo- { '> 

In Eq. (17), the Bose field obeys periodic boundary conditions, i.e., b((5) = b(0), and the Fermi 
fields obey the following boundary conditions: 

ai(/3) = ioi(Q) , a *(/3) = -ia*(0) , 
Pi(0) = i 0i(O) , ft(0) = -ift(O). 
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The free action for the single mode bosonic field Ssoib) is given by 

S B o(b) = f^dr b*(r) (d T - u ) 6(r) . 



(19) 



Then we can write the action S of the full fermion Dicke model, given by Eq. (11), using the free 
action for the single mode bosonic field Ssoib) defined by Eq. (19), plus an additional term that 
can be expressed in matrix form. Therefore the total action S can be written as 



S = S B0 (b) + I drY, p\{r) M(b*,b) p t (r) , 
the column matrix Pi(r) is given in terms of Fermi field operators in the following way 



(20) 



Piij) = 



OH T 



(21) 



and the matrix M(b*, b) is given by 

/ 

M(b*,b) = 



-1/2 



V 



(N)-^( 9l b(r)+g 2 b* (r) 



9i b* (r) + g 2 b (r) 



(22) 



/ 



the operators L and are defined by d T + 0/2 and d T — fl/2 respectively. Substituting the action 
5 given by Eq. (20) in the functional integral form of the partition function given by Eq. (17) 
we see that this functional integral is Gaussian in the Fermi fields. Now, let us begin integrating 
with respect these Fermi fields, therefore we obtain 



Z = J[dr]{b)}e SB0 (det M(b*,b)) N , 



(23) 



in this case, [dr](b)] is the functional measure only for the bosonic field. With the help of the 
following property for matrices with operator components 



det 



A B 
C D 



j = det (AD - ACA^B) , 



and determinant properties, we have that 

det M (&*, b) = det (LL*) det (l - N^L; 1 b + g 2 b*)^ 1 b* + g 2 b)) 
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(24) 



(25) 



Substituting Eq. (23) and Eq. (25) in Eq. (17), we have that 



Z f[dr)(b)} e?«> ' (26) 

with Za defined by 

Z A = J [d V (b)] exp (s B0 + Ntr\n(l- N~ l L~ l ( 9l b + g 2 b*)^ 1 ( 9l b* + g 2 . (27) 

We are interested in knowing the asymptotic behaviour of the quotient -J^- in the thermody- 
namic limit, i. e., N — > oo. With this intention, we analyse the asymptotic behaviour of the last 
defined expression Z A - First, let us scale the bosonic field by b — > y/Nb and b* — > y/Nb*, so that 
we get 

Z A = A(N) J[d7i(b)]exp(N$(b*,b)), (28) 

with the function $(6*, b) defined by 

$(&*, b) = S B0 + tr ln(l - L- 1 ( 9l b + g 2 b*)^ 1 ( 9l b* + g 2 ft)) . (29) 

The term A(N) in Eq. (28) comes from transforming the functional measure [dr](b)} under scaling 
the bosonic field by b — > \J~N b and b* — > \fNb*. The asymptotic behaviour of the integral 
functional appearing in Eq. (28) when iV — > oo, can be obtained by using the method of steepest 
descent [24]. In this method, we expand the function $(6*, b) around the point 6(r) = &o( r ) an d 
b*(r) = &o(t), which can be of two kinds. One kind that makes Re(Q(b* ,b)) maximum, and 
the other kind is defined as saddle point. We consider the first terms of the expansion in the 
integral functional, which are the leading terms for the value of the integral function. We can 
find the maximum points, or saddle points, finding the stationary points. The stationary points 
are solution of the following equations - ffey^ = and S f^*^ = 0. For the full Dicke model, the 

stationary points are constant functions b(r) = b and b*(r) = 6g. It is not difficult to show that 
for (3 < (3 C the stationary point is given by b = b^ = 0, which is a maximum point. The critical 
value (3 C is obtained by solving the following equation 

"° n tanhfMV (30) 



(gi + 92 



,2 



In this last equation, it is possible to find some solution for /3 C , in the case of (gi + g 2 ) 2 > uqQ. 
With this condition the system undergoes a phase transition. When the system has j3 < (3 C we say 
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that the system is in the normal phase. For (3 > (3 C the stationary points b{r) = b and b*(r) = 6g 
satisfy the following equation 



Phase transition happens if it is possible to find some real solution for |6 | 7^ in Eq. (31). It 
is only possible when (g\ + g 2 ) 2 > Uq^I and > j3 c . In the case of g\ 7^ and g 2 = 0, and 
also in the case of g\ — and g 2 7^ 0, the maximum points are a continuous set of values given 
by the expression b = pe 1 ^ and 6q = pe^ %<t> with e [0, 2ir) and p = |6 |, with |6 | defined by 
Eq. (31). In the case of g\ 7^ and g 2 7^ 0, we have two maximum points, which are given by 
^0 = ^0 — =•= j fc»o 1 5 with |6 | defined by Eq. (31). When the system has j3 > (3 C we say that the 
system is in the superradiant phase. 

Let us continue, with the computation of the asymptotic behaviour for the integral functional 
appearing in Eq. (28), for the thermodynamic limit, N — > 00. In following steps, we shall find this 
asymptotic behaviour when we only have one maximum point defined by b = b^. The resulting 
expressions will be useful for the normal phase of the full Dicke model, and also for the superradiant 
phase in the case of g\ 7^ and g 2 7^ 0. We consider the two first leading terms in the integral 
function appearing in Eq. (28) coming from the expansion of $(6*, b) around the maximal value 
6q = b , this expansion is given by 




(31) 



with defined by 




(32) 



$(&*, 6) = $(&*, 6 ) + 




dn dr 2 (6*(n) - b* , Kn) - 60) M$ 




the matriz M$, is given by 




(34) 



Substituting this expansion given by Eq. (33) in Eq. (28) we obtain 



ry 

Z A = e 




( 



b*(r 2 ) 
b(r 2 ) 



) 



) 



(35) 
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to obtain the last expression, we have applied the transformation b(r) — >\b(T) + boj/vN and 
b*(r) ^(b*(r) + 

K ) / m the functional integral involved. In order to make easier the integration 
of the functional integral given by Eq. (28), let us use the following transformation 

c(r) = a(g 2 b(T)+ gi b*(T)) 

c*(t) = a( gi b(r)+g 2 b*(r)), (36) 

the parameter a defined by the equation a 2 = (g 2 — <7i) -1 - It is worth mentioning that, the 
Jacobian of this transformation is 1. Applying this transformation in Eq. (28) we obtain that 

Z A = A(N) |[d77(c)]exp(7V$ 7 (c*,c)), (37) 
the function $j(c*, c) is given by 

$/(c*,c) = a 2 J drfacfr) - g 2 c*(r)) x 

x (d T - u> ) (<?! c*(r) - g 2 c(r)) + tr ln(l - a' 2 L; 1 c* L" 1 c^j . (38) 

The maximum point corresponds to Cq = c = ot{gi + g 2 )bo, the point 6q = b corresponds to a 
maximum for the function Re($(b*, b)). Using the same expansion given in Eq. (33) for $f(c*,c) 
and substituting in Eq. (37) we obtain that 

Z A = e ™(W J [dr]{c)] exp ^ £ dTl dT2 ( c *( Tl) > c(ri) ) M$J ^ ^ j j , (39) 
we have used the identity $j(cq, c ) = $(&o, M> an d the matrix is defined by 



(5 2 3>j(c*,c) gjj(cV) 

Jl/f _ I <5c*(ri) 5c*(r2) <5c*(ri) <5c(t2) 
A/l1 '' ~ 1 <5 2 $ 7 (c*,c) 5 2 $ 7 (c*,c) 

<5c(ri) <5c*(t2) 5c(ti)5c(t2) 



(40) 

C*=C=CQ 



At this level, it is convenient to use Fourier representation of the field c(r) in the functional 
integral Eq. (39). From boundaries conditions of the bosonic field b(r) and from Eq. (36), we 
deduce that c(r) and c*(r) satisfies periodic boundary conditions c(/3) = c(0) and c*(/3) = c*(0) 
respectively. Therefore Fourier representation of c(r) and c*(r) are given by 



c(r) = ^EcHe 4 -, 

c*(r) = 4lE^He^, (41) 
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the parameter oj takes the values: 2irn/f3, with n being all the integers. These values correspond 
to the Matsubara frenquencies for bosonic fields. Substituting this Fourier representation, Eq. 
(41), in Eq. (39), we obtain that 



Za = e N*(b* M J[d v ( c )}exp 
with <5 2 $(co>i, w 2 ) being defined by 



C (Ui) , C[U)i 



^2 j 



c*(to 2 ) 
c(w 2 ) 



(42) 



5 2 $n(wi,w 2 ) 5 2 <$> 12 (u u u 2 ) 
5 2 $ 2 i(wi,w 2 ) <5 2 $ 22 (wi,w 2 ) 

and each component of this matrix satisfies 



(43) 



5 2 $n(wi,w 2 ) 
5 2 $i 2 (wi,o; 2 ) 
5 2 $ 2 i(wi,o; 2 ) 
(5 2 $ 22 (wi, w 2 ) 



i • 5 2 $t(c* c 



1 r/ 3 

- / dndr 2 iri 
p Jo 



5c*(t 1 )5c*(t 2 
5 2 ®i(c*,c) 



-IOJ2T2 



C=C'=CQ 



1 

fiJo 







dTidr 2 e 1 ' 



5c* (ti) 5c(r 2 
5 2 $/(c*,c) 



l /■£ 

JO 



5c(ri) <5c* (r 2 
5 2 $/(c*,c) 



c=c'=cq 



C=C*=CQ 



/3 



5c(ri) &(r 2 ) 



(44) 



C=C*=CQ 



In this Fourier representation of the functional integral given by Eq. (42), the integral measure 
[c?77(c)] takes the tractable form l\ujdc*(u) dc*(u>). Using the expression for <3>/(c*,c) given in Eq. 
(38), we can calculate the matriz 5 2 <&(ui, u 2 ) with components given by Eq. (44). Performing 
these calculations we obtain that 



<5 2 $n(c<;i, u 2 ) = <5 2 $i 2 (wi,w 2 ) = S uij - U2 R(ui) , 

«J 2 $2i(wi,w 2 ) = 5 2 $ 22 (wi,w 2 ) = 6 ultUa S(u 1 ), (45) 
with S Ul jW2 being the delta Kronecker and the functions R(u>) and 5 (a;) are given by 



R(u) 
S(u) 



2 w fi'i #2 a 



iw 1 - 



Q 2 -Q 2 )a- 2 , //3fi A ~ 
tanh 



2fi A (w 2 + fi 2 i ; 



n A (u 2 + q a ; 



tanh 



2 
+ 



" Wo(ft + j2)ft+ vTfiij tan HT 



(46) 
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The expression for f^A is given by Eq. (32). Substituting the matriz 5 2 §(ui, 002), with components 
given by Eq. (45), in the functional integral appearing in Z A , given by Eq. (42), we obtain that 

Z A = e N ^ b oM J [ d?7 ( c )] exp ^ (^ S (uj) c(u>) c*(uu) + l - R(u>) (c(w) c(-uu) + c*(u>) c*(-u>) )) . (47) 

Performing this Gaussian functional integral, we finally obtain that 

z = e N*(b* ,b ) 2 _ll TT (2nif 

(S 2 (0) - R 2 (0)) 1/2 11 S (u>) S(-u>) - R 2 {u) • 1 } 

In order to find the asymptotic behaviour of when N — )■ 00, we must calculate / [drj(b)} e Sm 
appearing in Eq. (26). Using the free bosonic action Sbo given by Eq. (19), we obtain that 

\dr)(b)} e SB » = TT 2n% . (49) 

Substituting Eq. (48) and Eq. (49) in Eq. (26) we have that 

A = e NHKM * TT 1 (50) 

Z (H(0))W 11 //(c) ' ^ 

in this last equation, the function H{u) is given by 

HM= S(.)S( 7 )-R^) (51) 

W + (jJq 

The Eq. (46) gives the expresions for the functions S(u) and R(co), substituting these functions 
in Eq. (51) we obtain that 

/ x (g 2 -g 2 ) 2 n 2 , 2 //3fi A \ 

fi A (w 2 + ^a) ( W + w o) V 2 / 

2(fl 2 -fl 2 )^ 2 - (gf + gl)(n 2 + nl)co + 2 gi g 2 (g| - ^ 2 ) c ,//^a\ 
+ fi A (c 2 + Sft) (c 2 + c 2 ) tanh { 2 J' {b2) 

The expression, given by Eq. (50), with H(u) given by Eq. (52), provides a valid expression for 
the quotient in the normal phase, and also in the superradiant phase for the particular case of 
gi ^ and g 2 ^ 0. 
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5 Normal phase: j3 < f3 c 

In the normal phase, /3 < f3 c , from Eq. (31) we have that bo = b$ = 0, i.e. = Q. Subtituting 
this equality in Eq. (50) and Eq. (52), we obtain that 

z _ i n 1 

where 

ifj(w) = 1 + — ^tanh — + 
(w 2 + fi 2 ) (w 2 + wg) V 2 / 

+ 2(g? -^:;-^ + f ""° tanh ) . ,54) 
(uj 2 + IL 2 ) (u 2 + wg) V 2 / 

Making the analytic continuation (iu; — >■ i?) in Hi(u), we solve the equation Hj(-iE) = 0, which 
corresponds to the collective spectrum equation. Solving the equation, we have that 

2E 2 = ul + n 2 + 2 (<? 2 - £ 2 2 ) tanh + 

± ^ - + 4 (g\ (a* + fi) 2 - g\ (a* - fi) 2 ) tanh (^)) ^ ■ (55) 

It is interesting to see, when (3 = (3 C we find the following roots [13] 

£?i = (56) 

and 

E ( g 1 (n + g 2 + g 2 (n-^) 2 \ 1/2 

2 v (01 + 02) ; 1 ' 

With Eq. (55) we can obtain the collective spectrum for the two following known cases: the first 
one, when g 2 = 0, corresponds to the Dicke model considering the rotating wave approximation 
[10]. Here we have that 

2E = tu + tt ± ^(w -^) 2 + 4# 2 tanh(^^ . (58) 
The second one corresponds to gi — gi — g. Here we have that 

2E 2 = ul + Q 2 ± (^u 2 -Q 2 Y + 16g 2 u ntanh(^yj . (59) 

In the case of quantum phase transition, we are in the particular case where (5 = oo. Here, the 
collective spectrum corresponds to the Eq. (59) where: tanh(/3f2/2) = 1 [20]. 
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6 Superradiant phase: j3 > f3 c 



6.1 Case of g\ ^ and ^ 0: 



In the superradiant phase /3 > /3 C , in the case of g\ ^ and g 2 ^ 0, we have two maximum points. 
Both maximum points contribute equally to the partition function. Therefore J^, given by Eq. 
(50) must be multiplied by a factor 2. In this case b ^ 0, i.e. tt A 7^ ^- From Eq. (50), we have 
that the expresion for is given by 



— = 2e N * - TT - 



(60) 



where the factor 0, is defined by 



= " 



Q 2 



4 (gi + £ 2 ) 5 



+ In 



cosh 



/3n £ 



cosh (<$) 



(61) 



The function H n (u) has the form 



H n {u) 



1 



x 



u 2 + n\) (w 2 + ui 

2 , o2 , 2 (^1-^1) 



X 



n | 2 , gf 2 2 / 2 
C0 il\ U) + ; — CUq (il A - it ) 



(gi + 92? J (gi + 92) 2 

and setting u> = in Eq. (62), we obtain the expression for Hu(0), so that 

a 9i92 (n 2 A -n 2 ) 



(62) 



H n (0) 



(g 1 + g2 yni 



(63) 



Making the analytic continuation (iu — > E) in H n (u) given by Eq. (62), we solve the equation 
H n (—iE) = 0. The set of solutions, E, are the collective spectrum in the superradiant phase for 
the case of g\ ^ and g 2 ^ 0. Therefore, solving the equation we have that 



2E 2 



(91 + 92)' 



± 



u 2 + Q 2 A + 2 



( 9 2 i ~ 9l 
( 9i + 92) 



16gi #2 

! #1 + #2) 



2 W 



K - ^ 2 ) 



1/2 



(64) 
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For the particular case of g\ — gi — g, the collective spectrum energy takes the particular form 



2 E 2 = co 2 + n 2 A ± ((co 2 - n 2 A ) 2 + 4 co 2 n 2 ) 1 . (65) 

In the limit of zero temperature, /3 — > oo, from Eq. (31) we have that Qa = 4g 2 /co>o- Consequently, 
at zero temperature we obtain that [20] 

1/2 

(66) 



IE 1 



w + 



16 g 4 

i , 2 



± 



COn 



16 g 4 



UJn 



+ 4w 2 n 2 



6.2 Case of gi ^ and g 2 = 0: 

Now let us study, the case of rotating wave approximation, i. e., the case of g\ ^ and gi = 0, 
in the superradiant phase. Here, the expression for is obtained setting g 2 = in Eq. (28) and 
Eq. (29), therefore we have that 

Z A = A(N) j[dr](b)}e W (N^ gi (b*,b)), (67) 

the function <& gi (b*,b) is defined by 

<Z> gi (b*,b) = ' dTb*(u)(d T -u )b(u)+trln(l- g\ L~ x b L' 1 b*^j . (68) 

In last equation, Eq. (68), we can see that the function is invariant by transformation 

6(t) — > exp (i 9 t) b{r) and b*(r) — > exp (—i 9 r) b* (r) , where 9 is an arbitrary factor independent 
of r. This continuous invariance is responsible for the appearing of Goldstone mode in the system. 
In order to perform the functional integral given by Eq. (67), let us separate the function 6(r) in 
the following form 

6(r) = b c + V(r), 

b*( T ) = K + b'*(r), (69) 

where b c is a constant function, and the fields b'(r) and b'*(r) satisfy the following boundaries 
conditions fo'(0) = b'((3) = and b'*(0) = b'*((3) = 0. Using the representation b c = pe 1 ^ and 
b* = pe~ l< ^ in the functional integral given by Eq. (67) and Eq. (68), and after applying the 
transformation b'(r) ->■ e^fe'(r) and V*(t) ->■ e _ ^6'*(r), we obtain that 

/"CO /■ 

Z A = 2vr^(iV) ^ rfp 2 J [^(6 / )]exp(iV$ ^?1 (p,6 , ^6 , )), (70) 
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the function $ Sl (p, b'*, b') is given by 



%Mb'*,b') = jf dr(p + b'*(r))(d T -uj )(p + b'(r)) + 



(71) 



In the integral function appearing in Eq. (70), one variable of integration is p 2 . Here we use the 
steepest descent method in order to analyse the limit N — > oo, we find the stationary point with 

respect to the variable p 2 . Therefore, the stationary point satisfies the following equation -Jfy 

°\p ) p =p 

with b'*(r) = b'(r) = 0. In this case the value for po is the same as bo defined by Eq. (31) setting 
#2 = 0. Let us consider the two first leading terms in the functional integral appearing in Eq. (70), 
coming from the expansion of b'*, b') around the point defined by p and b'*(r) = b'(r) = 0, 

giving the maximum for Re($ gi (p, b'*, &')). This expansion is given by 



&'*,&') = <!>,(/,„. 0.0) 



p=p ,b'=b'*=0 



+ 



2 5(p 2 ) 2 



(p 2 -pI) 2 + 




(72) 



here the matriz M$ n is given by 

5 2 * 



91 

<56'*(ti)<56'*(t 2 ) 
5b' (n) Sb'* (r 2 ) 8b' (n) Sb' (t 2 ) 



56'*(n)56'(r 2 ) 



p = Po ,b'=b'*=o 



(73) 



Using this expansion given by Eq. (73) to perform functional integral given by Eq. (70), we have 
that 



Z A = 2vriv / iVe A ^ / 



VnpI 





dy e 



1 * 2 *31 
2 



p=p ,b'=b" 



X 



X 



J [drj(b')} exp 1 jf dn rfr 2 (^(n) , ft'fa)) M, 




(74) 



The expression 9l corresponds to the expression of defined in Eq. (61) taking g 2 = 0. The factor 
\/~N appearing in Eq. (74) come from the scaling p 2 — > p 2 / \fN . For iV — > oo, integrals appearing 
in Eq. (74) are Gaussians. We represent the functions b'(r) and b'*(r) in Fourier series, which do 
not possess the zero mode, since they satisfy the boundary conditions given by b'(0) = b'((3) = 
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and b'*(0) = b'*(f3) = 0. Therefore, performing the funtional integral and substituting in Eq. (26), 
we obtain that 

| = v ^ e «»^ n 1 , (75) 

functions (p 9l and Hu(lu) are given respectively by Eq. (61) and Eq. (62) setting g 2 = 0, and A 
is given by 

A _ 9i L \ h ( ) 

Making the analytic continuation (iuu — > E) in H n {u) given by Eq. (62) setting g 2 = 0, the 
collective spectrum is obtained by solving the equation H n (—i E) = 0. So that, we obtain the 
following spectrum 

£i = 0, (77) 

and 

El = col + n l + 2 w • (78) 

The particular value of the spectrum given by i?i = in Eq. (77) corresponds to the Goldstone 
mode [10]. 



6.3 Case of g\ = and gi 7^ 0: 

Now let us study the case of g\ — and g 2 7^ 0, in the superradiant phase. Here, the expression 
for is obtained setting gx = in Eq. (28) and Eq. (29). For this case we have 

Z A = A(N) J[dr](b)}e W (N<S> g2 (b*,b)), (79) 
the function <3> 92 (6*,6) is defined by 

^ g2 (b*,b) = jf" drb*(u) (d T -u )b(u)+tr\n(l - g\ L~ x b* L _1 bj . (80) 

In last equation, Eq. (80), we can see that the function Q g2 (b*,b) is invariant by transformation 
6(t) — > exp (i 6 r) 6(r) and b*(r) — > exp {—i 9 r) b* (r) , where is an arbitrary factor independent 
of t. This continuous invariance is responsible for the appearing of Goldstone mode in the system. 
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Since Eq. (80) is very similar to Eq. (68), we can see that, the calculation to obtain -J^ in the case 
of gi = follows the same steps as the calculation performed to obtain in the case of rotating 
wave approximation. Consequently, we have that 

where H n (uu) is given by Eq. (62) setting g 1 = 0, the value (fi g2 corresponds to the expression for 
4> defined in Eq. (61) taking gi = 0, and A is given by 

* (l- 4%- 'V . (82) 



Ha V 71 " P u o V sinh(/?r2A) 

Making the analytic continuation (iu> — > E) in H n (ou) given by Eq. (62) setting gi = 0, the 
collective spectrum is obtained by solving the equation H n (—iE) = 0. So that, we obtain the 
following spectrum 

£i = 0, (83) 

and 

El = u 2 Q + Q 2 A - 2 cuo n . (84) 

The particular value of the spectrum given by E 1 = in Eq. (83) corresponds to the Goldstone 
mode. 



7 Summary 

In this paper, using the path integral approach and functional methods, in the thermodynamic 
limit N — > oo, we find the asymptotic behaviour of the partition function and collective spectrum 
of the full Dicke model in the normal and superradiant phase. In our study we distinguish three 
particular cases. The first one corresponds to the case of rotating wave approximation, g 1 ^ and 
gi = 0, in this case the model has a continuous symmetry, which is associated to the conservation 
of the sum of the number excitation of the N atoms with the number excitation of the boson field. 
The second case corresponds to the model with g 1 — and g 2 ^ 0, in this case the model has a 
continuous symmetry, which is associated to the conservation of the difference between the number 
excitation of the iV atoms and the number excitation of the boson field. The last one, corresponds 
to the case of g\ ^ and g 2 ^ 0, which the model has a discrete symmetry. The phase transition 
in each case is related to the spontaneous breaking of their respective symmetry. In the case of 
rotating wave approximation, and also in the case of g 1 = and g 2 ^ 0, the collective spectrum has 
a zero energy value, corresponding to the Goldstone mode associated to the continuous symmetry 
breaking for these cases. 
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